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ABSTRACT 


The  siruLlural  response  of  three  full -scale  helicopter  rotors 
has  been  examincrl.  Airloads  predicted  by  CAMRAD  were 
compared  by  htirmonics  with  test  measurements.  Rapwise 
bending  moments  were  calculated  from  the  CAMRAD  air¬ 
loads  using  three  methods;  force  integration,  curvature,  and 
a  fmite-diffcrcnce  method.  The  force-integration  and  cur¬ 
vature  moments  were  calculated  with  the  CAMRAD  pro¬ 
gram  itself;  however,  the  finite-difference  metliod  was  cal¬ 
culated  externally  by  a  forced  response  program,  Blade  Rc- 
sjxmse  to  Aeiudynamic  Loading  (BRAL).  1  nc  BRAL  anal- 
>  sis  w  ith  measured  airloads  was  shown  to  agree  well  with 
sdain-gaugc  measurements  on  the  CH-34  rotor  in  a  wind 
tunnel.  The  CAMRAD/force-integralion  bending  moments 
were  shown  to  accumulate  error  as  the  integration  pro¬ 
gressed  from  lip  to  root.  The  CAMRAD/curvature  mo¬ 
ments  agreed  well  with  the  finite-difference  moments  over 
most  of  die  blade;  however,  at  regions  of  rapidly  chang¬ 
ing  stiffness  the  agreement  was  poor.  When  applied  to  the 
BRAL  solution,  force  integration  was  shown  to  give  excel¬ 
lent  results,  provided  the  integration  had  a  .small  step  si/c. 
Shear  forces  calculated  by  CAMRAD  and  BRAL  showed 
fair  agreement. 


INTRODUCTION 

fhe  calculation  of  helicopter  rotor-blade  bending  momenus 
IS  acroelastic  in  nature.  Aerodynamic,  inertial,  and  clas¬ 
tic  forces  must  balance  to  provide  equilibrium  of  the  ro¬ 
tor  blade.  .Although  the  Ixmding  moments  are  embedded 
m  ’.he  cquilmriuin  solution,  they  are  typically  calculated 
separately  after  die  equilibrium  solution  has  been  obtained. 

I  here  are  two  common  ways  of  doing  this.  The  bending 
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moments  are  cither  calculated  as  the  integral  of  the  aerody¬ 
namic  and  inertial  forces  or  diey  are  calculated  as  the  prod¬ 
uct  of  stiffness  and  curvature. 

The  relative  iiieriLs  ol  diese  two  methods  have  been  exam¬ 
ined  in  the  literature.'"’'  All  have  eoiieluded  diat  force  in¬ 
tegration  is  more  accurate  and  that  it  converges  with  fewer 
modes.  Biclawa^  pointed  out  that  the  incrca.sed  accurai_\ 
of  force  integration  is  paid  for  by  increased  complexity 
and  diat  force  integration  requires  greater  attention  to  de¬ 
tails  of  its  implementation.  Numerical  problems  arise  for 
both  methods.  The  curvature  method  requires  the  calcu¬ 
lation  of  the  second  derivative  of  die  dellected  shajie.  aiul 
the  force-integration  method  requires  that  the  difference  of 
large  forces  be  integrated  along  the  blade  span.  If  the  struc¬ 
ture  has  rapid  changes  in  stiffness  and  mass,  there  must 
also  be  rapid  changes  in  the  dellected  shape.  This  requires 
many  modes  with  a  mtxial  analysis  or  many  elements  w  ith 
a  (initc-elemeni  analysis.  This  also  requires  that  the  force- 
integration  niediod  have  small  enough  intervals  to  accu¬ 
rately  evaluate  the  della'tions  and  die  spanwise  distribution 
of  structural  properties. 

Hodges’'  'leveloped  the  direct  method  of  Ritx,  as  a  liiiite- 
clemcnt  method  and  demonstrated  that  it  could  achieve  die 
‘'exact"  solution  w  hen  elements  spanned  regions  of  analytic 
stillness  and  mass.  He  compared  diis  with  conventional  ap¬ 
plications  of  the  Rayleigh-Rit/.  niediod  in  which  both  force 
integration  and  curvature  were  used  to  calculate  Ixmding 
moments.  He  concluded  that  die  force-integration  mediod 
was  more  accurate  than  the  curvature  mediod  but  that  die 
Riu  finite-element  method  was  superior  to  both. 

In  correlating  results  from  a  Comprehensive  Analytical 
Model  ol  Rotorcraft  Aerodynamics  and  Dynamics.  CAM¬ 
RAD,  with  llight  test  data  for  the  ,SA.14‘V2  rotor  it  was 
reported  that  die  calculated  llapwise  bending  monieiit.s 
were  .sensitive  to  die  blade-mass  distribution.'’  At  that  lime 
CAMRAD  calculated  bending  moments  using  the  force- 
integration  metluHl.  CAMRAD  was  changerl.  hovvever.  to 
calculate  bending  momenis  using  the  curvature  method;  the 
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result  was  a  signilicant  improvement  in  loads  correlation.^ 
Subsequent  loads-correlalion  elTorLs  using  CAMRAD  have 
relied  on  the  curvature  method.^'’  The  superiority  of  the 
curvature  methtxl  in  CAMRAD  is  surprising  considering 
the  conclusions  of  Rel’s.  1  -4. 

This  paper  examines  the  calculation  of  llapwisc  bending 
moments  alter  an  equilibrium  solution  has  been  obtained. 
CAMRAD  is  used  to  obtain  the  aeroelastic  equilibrium  so¬ 
lution  lor  three  full-scale  rotors:  the  CH-34,  SA349/2,  and 
BO- 10.3.  Subsequent  llapwisc  bending  moment  calcula¬ 
tions  will  be  made  with  CAMRAD  using  both  the  force- 
integration  and  curvature  mcthmls.  The  airloads  from  the 
CAMRAD  equilibrium  solution  arc  also  used  as  the  forcing 
function  for  an  independent  fmite-diffcrcncc  response  anal¬ 
ysis  dial  is  used  to  evaluate  results  of  the  two  metliods  cal¬ 
culated  in  CAMRAD.  The  finite-difference  solution  is  used 
to  make  general  conclusions  about  tlic  force-integration  and 
curvature  methods  for  calculating  bending  moments  for  re¬ 
alistic  full-scale  helicopter  rotors. 


CAMRAD  CALCULATION 

CAMRAD  is  used  to  calculate  trim,  airloads,  blade  dellec- 
iions,  and  structural  loads.  A  detailed  description  of  the 
theoretical  basis  of  CAMRAD  is  given  in  Ref.  8.  The 
blatle  structural  analysis  follows  the  work  of  Houbolt  and 
Brooks'*  w  ith  important  nonlinear  effects  included  in  the  in¬ 
ertial  diid  aerodynamic  forces.  The  aerodynamic  model  is 
based  on  lifiing-linc  theory,  using  steady  two-dimensional 
airfoil  chiiractcristics  and  a  vortex  wake.  Blade  aerody¬ 
namic  and  structural  models  provide  the  section  foiccs. 
Hqudibriiim  of  inertial,  aerodynamic,  and  clastic  forces  dc- 
lermines  tlic  differential  equations  of  motion.  The  blade 
degrees  of  freedom  are  llap/lag  bending,  rigid  pitch,  and 
elastic  torsion  modes. 

Orcc  CAMR.AD  has  calculated  the  trim,  airloads,  and  blade 
dellections,  the  bending  moments  are  calculated.  The  orig¬ 
inal  version  of  CA  MRAD  dcKumcnted  in  Ref.  8  calculates 
the  bending  moments  using  the  force-integration  methtxJ. 
As  mentioned  m  the  introduction,  a  mrxiilication  was  writ¬ 
ten  to  CAMRAD  which  calculates  bending  moments  using 
the  curvature  methixl.  Results  from  both  methods  will  be 
shown. 

.Additional  modilicalion  to  Itie  CAMRAD  analysis  used 
here  has  iiiqiruved  the  high-speed  llighl  airload 
calcuhition.  "  An  improved  wtike  model  has  been  incor¬ 
porated  III  die  analysis.  It  includes  a  du;il-peak  circula¬ 
tion  model,  to  more  accurately  model  tlic  negative  loading 


at  die  tip  on  die  advancing  side  of  die  disk,  and  second- 
order  iifting-line  theory,  to  better  nuxici  blade-vonex  inter¬ 
action.  This  is  essentially  die  aerodynamic  model  available 
inCAMRAD/JA,'' 


FINITE-DIFFERENCE  STRUC¬ 
TURAL  RESPONSE  ANALYSIS 

To  aid  in  die  comparison  of  the  two  bending-moment  calcu¬ 
lations  implemented  in  CAMR.AD  a  line-resoluiion,  linite- 
differcnce  structural  response  analysis  has  also  Ix'cn  used. 
The  Blade  Resixinse  to  Aerodynamic  Loading  (BRAL)  pro¬ 
gram  originated  from  die  work  of  Esculier  and  Bousman.' ' 
BRAL  calculates  the  coupled  llap/lag/torsion  structural  re- 
spon.se  of  die  helicopter  rotor  to  given  aerodynamic  load¬ 
ing  as  described  in  Ref.  12.  A  periodic  solution  is  as¬ 
sumed  in  me  time  domain  and  is  represented  by  t-ourier 
harmonic  functions.  The  linear  Houbolt  and  Brooks'*  equa¬ 
tions  of  motion  arc  solved  for  each  harmonic  as  a  two- point 
boundary  value  problem.  In  effect,  this  is  a  linite-difference 
methtxl  with  several  hundred  span  wise  grid  poiiiLs,  The 
BRAL  analysis  has  been  validated  using  airloads  measured 
on  the  CH-34  in  a  wind  tunnel.'^ 

These  measured  dirtoetd.t  wore  used  as  the  forcing  func¬ 
tions  to  calculate  bending  moments  using  BR.AL  with  400 
spanwisc  intervals.  Figure  1  compares  the  BRAL  llap- 
bending  moment  with  wind-ninncl  measurements.  This  fig¬ 
ure  shows  the  flap-bending  moment  time-history  at  45Cf 
blade  span  in  the  upper  left-hand  corner.  Following  the 
lime-history  arc  spanwisc  distributions  of  the  steady  and 
first  six  harmonics  of  the  llap-bcnding  moment.  The  am¬ 
plitude  scale  of  the  harmonic  components  was  allowed  to 
vary  so  dial  the  spanwisc  shape  of  the  components  could  be 
seen.  Although  BRAL  uses  a  linear  structural  model,  the 
resulting  correlation  is  excellent. 

In-vacuum  blade  modes  were  calculated  to  ensure  diat  die 
structural  models  were  as  simiLir  as  possible.  CAMR.AD 
modes  were  calculated  using  the  llutlcr  analysis  option, 
and  BRAL  modes  were  calculated  by  a  companion  pro¬ 
gram,  Coupled  MtxJcs  (CMODE.S).  CMODE.S  calculates 
the  rtHils  to  the  homogcncou:;  form  of  die  equations  of  mo¬ 
tion  found  in  BRAL,  using  the  frequency  determinant." 
Frequency  diagrams  for  die  CH-34,  .SA349/2,  and  130-10.3 
rotors  are  shown  as  Figs.  2,  3,  and  4,  resjx'ctively.  Agree¬ 
ment  is  excellent  for  die  CH-34  rotor  and  fair  for  the  otli 
crs.  The  greatest  disagreement  ixcurs  for  the  first  torsion 
frequency.  The  CMODE.S  program  dtx's  not  include  steaily 
deformation  caused  by  center  ol  gravity  aiui  tension-center 
offsets  from  the  ehistic  axis.  When  these  offsets  tire  set  to 


zero  or  when  zero  bending  degrees  ol  freedom  are  used  for 
the  menn  delleeted  shape  in  CAMRAD,  the  torsion  frequen¬ 
cies  show'  much  better  correlation. 


ROTOR  AIRLOAD 
CALCULATIONS 


As  described  earlier,  tlie  calculation  of  bending  moments 
begins  with  the  acroelastic  equilibrium  solution.  Normal 
airloads  calculated  by  the  CAMRAD  equilibrium  solution 
are  presented  in  this  section  with  test  measurements,  where 
available.  It  is  important  to  note  tliat  the  CAMRAD  air¬ 
loads  arc  identical  for  each  of  the  BRAL,  CAMRAD/forcc- 
integration  and  CAMRAD/curvaturc-monicnl  results.  Al¬ 
though  not  shown,  chordwise  airloads  and  pitching  mo¬ 
ments  die  coiiiiiion  to  tin-sc  stru..mral  moment  calculations 
as  well. 

Figures  5-7  show  the  blade  nonnal  airloads  calculated  by 
CA.MRAD  and  measured  in  tests,  with  the  c.xccptton  of  the 
BO- 105,  for  which  no  test  data  arc  available.  The  CH-34 
and  S  A349/2  rotors  arc  shown  at  a  similar  high-speed  flight 
condition;  the  BO-l()5  is  at  a  lower  speed.  Wind-tunnel 
test  measurements  of  the  CH-34  agree  very  well  with  the 
steady  component  of  CAMRAD  airloads.  Agreement  is 
fair  for  the  first  and  second  harmonic  components,  but  har¬ 
monics  components  alxfvc  the  second  harmonic  show  less 
agreement.  The  correlation  between  the  SA349/2  flight- 
test  measurements  and  CAMRAD  airloads  is  similar  to  the 
CH-34  results,  although  with  only  three  spanwise  measure¬ 
ments,  conclusions  arc  limited.  A  distinctive  similarity  may 
be  seen  in  the  CAMRAD  predicted  normal  airloads  for  all 
three  rotors. 


DEFLECTED  SHAPES 


The  equivalence  of  the  structural  models  of  BRAL  and 
CAMRAD  may  be  seen  (Figs.  8-10)  by  comparing  dc- 
(lected  shapes.  CAMRAD  has  iterated  to  calculate  the 
acroclastically  U'immcd  solution,  and  BRAL  has  calculated 
the  sU'uciural  response  to  the  CAMRAD  airloads.  The  de¬ 
flected  shapes  would  be  the  same  for  identical  structural 
models,  since  the  airloads  arc  identical.  Therefore,  differ¬ 
ences  seen  here  reflect  differences  in  the  structural  models. 


shown.  The  results  are  very  similar  for  the  three  rotors.  Dis¬ 
agreement  between  the  two  analyses  is  greatest  at  tlie  first, 
third,  and  fifth  harmonics,  which  arc  close  to  the  first,  sec¬ 
ond,  and  third  flap  modes.  This  is  particularly  noticeable  for 
the  first  htu'moiiic  cosine  component  ol  tfie  flap  deflection 
for  the  SA349/2.  The  pio.ximity  of  the  flap  modes  to  the  ro¬ 
tor  harmonics  is  determined  by  the  su-uctural  model.  It  was 
seen  in  Figs.  2-4  that  the  BRAL  (CMODES)  and  CAM¬ 
RAD  modes  arc  not  identical,  'fhercforc,  the  disagreement 
in  dcflcctcrl  shapes  may  be  due  to  small  differences  in  the 
modal  frequencies. 


BENDING  MOMENTS 


Figures  11-13  show  the  resulting  llapwise  bending  mo¬ 
ments.  Test  measuremeuLs  are  shown  with  BRAL.  C,\M- 
RAD/curvaturc,  and  CAMRAD/  force-integration  calcu¬ 
lations.  The  time-history  plot  of  Fig.  11  shows  C.\M- 
RAD/curvaiiirc  and  BRAL  piedictions  agreeing  well  with 
one  anotlier.  The  time-history  plots  of  Figs.  12  and  13 
show  fair  agreement  between  the  flight-test  measurements. 
BRAL,  and  CAMRAD/curvature.  The  CAMRAD/force- 
integration  results  show  poor  correlation  for  all  thicc  rotors. 

Examination  of  the  luirmonic  compvjiieius  of 
Figs.  1 1-1.3  reveals  two  consistent  discrepancies.  First,  the 
CAMRAD/forcc-integration  has  the  correct  zero-moment 
boundary  condition  at  tlie  tip  of  the  blade,  but  it  appears  to 
accumulate  error  as  forces  are  integrated  toward  the  root. 
This  is  most  obvious  for  the  zero-moment  boundary  con 
dition  at  the  flap  hinge  of  the  Cl  I  33  rotor  (Fig.  11).  This 
would  also  be  true  for  the  Sz\ 349/2,  but  the  BRAL  calcula¬ 
tion  stops  at  flic  outboard  lag  hinge,  enforcing  the  inboard 
flap  hinge  boundary  condition  assuming  rigid  txxly  motion. 


The  second  discrepancy  seen  in  Figs.  11-13  is  flie  hiige 
changes  bending  moment  predicted  by  tlie  C.'XM- 
RAD/curvalurc  mcfliod  at  the  root  end  of  the  blades.  This 
is  caused  by  large  changes  in  bending  stiflTess  ihat  occur 
at  the  root  of  all  three  blades  (Fig.  14).  Seven  llapwise 
bending  modes,  as  used  here,  aie  not  enough  to  capture  the 
rapid  changes  in  curvature  required  to  give  a  smooth  prod- 
iK'i  nf  curvaiiirc  and  s’ilTncys.  Tb"  phenomenon  a'so  ap 
pears  near  the  tip  of  the  Sz\  .349/2  rotor  blade  (best  illustrated 
by  the  first  harmonic  plot.  Fig.  12).  In  general,  it  is  seen  for 
all  three  rotors  that  the  results  of  the  CAMRAD/curvature 
mcthtxl  agree  well  wifli  those  of  BRAl,  away  from  rapid 
changes  in  bending  stiffness.  Ct/ 
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Figures  8-iO  show  the  llapwise  deflected  shapes  calculated 
by  BRAL  and  CAMRAD  for  the  CA.MRAD  airloads  just 


DISCISSION 

l  Niiii;  iJcDik'ai  airlojils,  (lilTi'rcikL's  li:i\c  Ivoii  shiiN.'-n 
111  liic  tx'nding  nioiiK'ni^  cak'ulalcd  by  iho  BKAL, 
CAM  RAD  Tore  i.'-iniL'graiion,  and  CAMRAD/eurvalurc 
laeiluKls.  To  die  degree  dial  die  dedeeied  shapes  are  siindar, 
diese  diri'erenecs  are  atiribuied  lo  inipleiiientalion  problems 
vutli  die  roree-integration  and  eurvature  mediods.  Depend¬ 
ing  on  ihe  radial  siadi'n  seleeied,  eidier  meduxi  will  show 
iiiaju:  errors.  Figure  1.x  deinonstraies  die  severiiy  of  die 
prohlein  near  the  blade  root  lor  both  meihtKis. 

Insight  can  be  gained  by  ealeulating  bending  moinents  di¬ 
rectly  from  the  BRAL  solution.  Thus  die  siaie-veetor  (de- 
ileeliori.  slope,  bending  moment,  and  shear)  solution  at  .M)() 
spanw  isc  intervals  was  used  to  caleulaie  the  ’  ending  mo- 
nienis  by  the  fcirce-iniegraiion  and  euiwatiire  methods.  Re- 
sulimg  bending  moments  may  he  eonipared  w  ith  die  bend¬ 
ing  morneni  of  die  state  vector.  The  slate  sector  was  cai- 
suLited  using  die  uncoupled  dap  equation  in  BR.AL,  for  an 
untwi-.ted  C11-T4  rotor  widi  the  CA.MR.AD  airloatis  shown 
in  f  ig.  .X.  Inertial  forces  were  ctilculaled  from  the  suite- 
vector  dellected  shape  assuming  simple  harmonic  motion', 
aerodynamic  forces  were  taken  from  the  C.AMR.AL)  solu¬ 
tion;  and  the  curvature  was  calculated  as  the  derivative  ‘>i' 
the  iiUUC-'vcvior  siojic.  the  spanwise  integration  for  the 
lorce-inlegraiion  method  was  calculated  using  the  tra[ie- 
/oitlal  rule. 

Figure  16  shows  die  moments  due  lo  the  indiv  idual  forces: 
aerodynaimc,  inertial,  and  centrifugal.  FTgure  17  shows  die 
IfR.AI.  state-vector  moment;  curvature-method  moment; 
lorce-inlegradoii  moment,  1(X)  intervals;  tiiul  the  force- 
n.tegralion  moment,  .T(X)  intervals.  .Mlhough  the  step  si/e 
variable  in  tne  BRAL  analysis,  foice-iniegraiion  nio- 
menlv  calculated  from  the  IfR.AL  stale  vector  aiul  m  C.A.M- 
R  \D  ii-e  uniformly  spaced  intervals.  ITie  results  seen  in 
big.  17  are  very  similar  lo  those  from  the  C-AMRAD/Toae- 
mtegration  and  CAMRAD/curvtitiire  methods.  The  force- 
iniegraiion  moments  with  100  intervtils  accumulates  error 
lo'.v  ard  the  blade  root,  and  the  curvature  method  is  in  error 
i:i  the  reeion  of  rapitlly  changmg  stillness. 

Ihe  lorce-inlegration  results  with  l(X)  tind  .slXi  intervals 
sho'vv  dial  an  insuflicieni  number  of  spanwise  integra¬ 
tion  steps  can  rcsml  m  di  '  kind  of  error  seen  in  the 
(  AMR.ADTorce-iniegralion  results.  This  iiulicaKa!  that 
a  r.', ill, lit  rotor  such  ;is  the  ('H-T4  rec|uires  at  least  tOO 
;  111  .'.  1  'tucrvals  using  this  iinplcmenltition  ol  the  lorce- 
m;!  method.  It  also  indicates  lluii  the  itiiegniiion 

1,  I  riii.al  whereas  numenc ''  piol'ieins  lated 

■Mill  111'.'  diiierc  iice  of  larue  miiiilx'rs  do  not  ctiiise  problems 
III  this  calciilat  loa. 


Figures  16  and  17  taken  logedier  illustrate  another  possible 
source  ol  error  for  the  lorce-iniegration  method.  Near  die 
blade  riKil,  the  magnitude  of  the  total  structural  inomeul  is 
less  than  IFf  olThe  kirgest  moment  component.  Therefore, 
a  small  error  m  die  calculation  of  aertxlynamic,  inertial,  or 
centrifugal  moments  may  give  rise  to  a  large  error  in  the 
smuctural  moment. 

.Additional  insight  may  he  gained  by  considering  the  indi¬ 
vidual  forces  calculated  with  die  force-integration  nielluxl. 
Figure  IS  shows  die  individual  component  foices.  and 
Fig.  Id  shows  the  sum  of  die  component  forces  compared 
with  the  BR.AL  state-vector  shear  force.  ITie  comporenis 
of  the  shear  lorce  and  die  sum  of  shear  forces  were  com¬ 
puted  at  100  intervals  along  the  blade.  .Agreement  between 
the  .sum  of  die  coiiiponeni  forces  and  the  siaie-v  ecior  shear 
force  is  very  good.  This  indicates  dial  fewer  intervals  are 
needex;!  to  accurately  calculate  shear  forces  than  are  needed 
for  bending  momeiiLs. 

1-igure  20  tienionstraies  die  agreement  Ix'iween  shear-lorce 
calculations  of  BR.AL  and  C'.AMRAD.  The  agreement  ol 
shear  forces  is  very  similar  to  the  tigreemcni  seen  m  the  de- 
llecied  shapes  in  Fig.  S.  This  is  not  surprising  since  the 
airloads  are  identical  for  the  two  analyses  and  the  inertial 
forces  aie  calctiiatcd  Irom  the  dede "led  shapes. 


CONCLUSIONS 

T  he  i.solatcd  proiTem  of  calculating  llapwise  bending  mo¬ 
ments  from  a  known  aerixtlasiic  equilibrium  solution  has 
been  c.xplored.  A  comprehensiv  e  roiorcrai't  analy  sis,  C.AM- 
RAD,  was  used  to  obtain  this  solutuin  and  lo  tlemoiisdaie 
the  force-inie;  ■  and  curvature  meiluxls  ol  calculat¬ 
ing  bending  nux  A  liniie-dillerence  method,  BR.AL, 
was  used  lo  obu..  .curtiie  llapwise  beiuling  moments. 
These  vvere  compared  with  calculations  using  ilie  force- 
iiitegraiioii  and  curvature  mcThexls  lo  make  general  ctinclu- 
sioiis  alxHil  die  two  methtxis. 

It  was  concitHlcd  that  the  C.AMRAD/forcc-integralion 
nielluxl  accumulates  error  as  the  iiuiiiieni  is  uiiegraled  Irom 
lip  lo  root.  The  C.AMRAD/curvaiure  meiiuKl  gives  results 
that  tire  gooil  or  bad  deivndmg  on  the  radial  station.  Raili.il 
stations  away  from  rajiid  changes  in  stillness  were  shown 
to  give  cxccilc.r  :  'sullv;  '■.■^i|i|,^  .ii.-ar  these  rapid  stillness 
changes  must  be  ignored.  TTie  C'.AMR.AD  shear  lorces. 
which  tire  the  integrated  aerixlynamic  tind  inertial  knees, 
showed  no  accumulalion  oi  ernw  :■■■■  mi  lU,-  hi-.a  ■ 

I  iiesc  rcstilLs  are  general  to  the  extent  that  they  are  caused 
by  the  complex  structure  of  existing  full-scale  roiors.  T  hat 


IS  noi  10  say  ihal  ihcsc  problems  cannot  be  overeome  by  an- 
uitiei  implemenuition  with  suHieient  degrees  ol  freedom  or 
by  a  more  sophisticated  integration  scheme.  In  fact,  when 
applierl  to  the  BRAL  solution,  the  foree-integratior;  r.iclhod 
was  slio'wn  to  give  accurate  bending-moment  results,  pro¬ 
vided  there  were  enough  integration  steps. 


APPENDIX:  DESCRIPTION  OF 
ROTORS  ANALYZED 

The  structural  resptrnses  of  three  rotor  systems,  the  CH-34, 
SA.349/2,  and  BO-105,  were  calculated  and  compared  with 
test  results.  The  source  of  the  structural  properties,  test  data, 
and  a  description  of  diniculties  encountered  in  modeling 
these  rotors  follows. 

The  structural  properties  for  the  CH-.^4  rotor  were  Liken 
from  Ref.  15.  Calculated  airloads  and  bending  moments 
were  compared  with  test  data  from  the  wind-tunnel  lest  de¬ 
scribed  in  Ref.  1 3.  Bending  moments  were  converted  from 
blade  stresses,*^  using  the  stress-moduli  value  reported  in 
R  'f.  15.  The  Cfl  .34  rotor,  with  coincident  llap  and  lag 
hinges  and  blade  designed  to  minimize  llap/lag/torsion  cou¬ 
plings.  was  easily  modeled  by  boili  CAMRAD  and  BRAL. 
The  wind-tunnel  test,  however,  had  a  split  tip-path  plane  be- 
caie  e  of  the  swashplate  actuator  insUillation.  This  was  ivn 
modeled  by  the  analyses,  which  assume  identical  blades. 

The  structural  properties,  airloads,  and  bending  moments 
for  the  SA349/2  rotor  were  all  uiken  from  liic  llight  test  re¬ 
ported  in  Ref.  16.  .Static  values  of  the  tlap-bcnding  mo¬ 
ment  were  not  used  in  comparisons,  a  result  cT  calibration 
problems.  The  CAMRAD  model  of  the  .S.-\349/2  hinge  ar¬ 
rangement  is  only  appro.xirnated.  The  hinge  sequence  is 
llap/feather/Iag  from  inboard  to  outboard;  therefore,  as  the 
hlade  feathers,  the  a.xi.s  of  flic  lag  hinge  changes  relative  to 
the  fixed  axis  of  the  lltip  hinge.  CAMRAD,  however,  re¬ 
quires  that  tfie  (lap  and  lag  hinge  axes  be  perpendicular  to 
one  another.  The  CA.\'RAD  and  BRAL  results  shown  here 
model  the  llap  hinge  accurately,  hut  the  lag  hinge  remains 
lixed  regardless  of  pilch  input. 

1  he  structural  properties  of  the  BO- 105  are  available  m 
.Refs.  17  and  lb.  Bending  moments  were  taken  from  the 
liight  test  of  Ref.  17.  Airf  /ad.s  have  not  been  measured; 
'he, -eh;.''''  '’f'  tnst  data  ..re  compared  with  CA.MRAD  pre¬ 
dictions.  Canlilcvci  ruit  end  Ixnmdary  conditions  are  used 
•  rnM/Cq  this  pimir'!.-'-': 


Rotor  parameters  of  interest  are  listed  in  Table  1.  Shown 
in  Fig.  14  are  tlic  llapwisc  stillness  and  mass  distributions 
for  these  rotors.  All  three  rotors  have  a  build  up  of  structure 
neru"  the  nxit  end  of  the  blade  seen  as  an  increase  in  stiffness 
and  mass.  Additionally  the  .SA349/2  has  several  changes  in 
stiffness  and  mass  near  die  blade  tip,  and  die  CH-.34  has 
small  changes  in  mass  distribution  owing  to  insirumenU'- 
lion  and  dp  weights. 
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Table  1  Rotor  parameters. 


Item 

CH-.T4 

S  A  349/2 

BO- 105 

Numlx;r  of  blades 

4 

3 

4 

Rotor  radius,  m 
(in.) 

8.53 

(336.00) 

5.25 

(206.85) 

4.91 

(193.32) 

Rotor  speed,  rjtm 
(rad/sec) 

221.68 

(23.214) 

387.00 

(40.526) 

425.(X) 

(44.506) 

Blade  chord,  m 
(m.) 

0.416 

(16.40) 

0.350 

(13.80) 

0.270 

(10.64) 

Solidity 

0.06220 

0.06366 

0.07010 

Fla))  hinge  offset,  epIK 

0.0357 1 

0.02095 

- 

Feathering  hinge  offset,  cv/R 

0.08(KK) 

0.0-1762 

0.03500 

l.ag  hinge  offset,  (;<;/R 

0.0357 1 

0.09048 

- 

Blade  cutout,  r/R 

0.1500 

0.2486 

O.IOOO 

Precone,  deg 

0.0 

0.0 

2.5 
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[  ig.  1  C'linipanson  ol  llapwisc  Ix-nding  momcnis.  CH  34  rotor:  =  0.39,  shall  angle  = 
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F  ig.  .5  Compon.son  of  normal  airloads,  CH-34  rotor:  =  0.39,  shaft  angle  =  -.S-’. 
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l  ig.  6  Coniiuinson  of  normal  airloads,  .SA349/2  rolor:  a  =  0.378,  shall  anolc  =  -5.: 


CAMRAD  AIRLOAD 


I  IMt  HISTORY  f  R  0  7b 


1ST  HARMONIC 


7ND  HARMONIC 


COSINfc 


SINE 


1000 


z 

n 

< 


bOO 


0 


5  ^00 


1000 


TRD  HARMONIC 


4TH  HARMONIC 


f)TH  HARMONIC 


COSINE 


SINE 


700 


100 


6TH  HARMONIC 


BKAL  CAMRAD  AiHLOADSi 
CAMRAU 

TIMf-  HISTORY  i  R  0  80 


IS  r  HARMOMIC 

cos:  NT  ,  SINt 


STEADY 


I  R 


2ND  HARMONIC 


I 


3RD  HARMONIC 


4TH  HARMONIC 


SIH  HARMONIC  BTH  HARMONIC 


I  T2.  S  {  onipariMiii  1)1  tlapwiso  dolloclL’d  shapes.  C'll-^4  rotor:  /j  =  O.dd,  shall  anelo  - 
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i  ii',  H  (  iiiiiiiorison  ol  lUipwisc  clctlcclctl  sliapcs,  SA349/2  rotor:  =  0.378,  shall  angle  = 
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1  lii.  10  (Omparison  of  Ikifiwisc  ticllctial  shapes,  HO- 103  rotor:  u  =  0.278,  shall  anple  = 
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Fig.  1 1  Comparison  of  llapwisc  bending  moinenls,  CH-.^4  rolor:  =  0,39,  shall  angle  =  -5°. 

•  FLIGHT  TEST 

-  BRAL  ICAMRAD  AIRLOADS! 

- CAMRAD  (CURV.) 

- CAMRAD  (F.l  ) 


TIME  HISTORY  r/R  0  46 


1000  L  -  -  -  1  ‘  1  i  -  L  _ 

0  45  90  135  180  225  270  315  360 

AZIMUTH,  cley 

1ST  HARMONIC 

COSINE  SINE 


STEADY 


0  2  4  6  8  1  0 

r/R 


i  lg.  1 2  ('omparison  of  flapwise  bending  moments,  .SA349/2  rolor:  /i  =  ().37(S,  shall  angle  =  -3.5°. 
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Fig.  14  Flapwisc  stiffness  and  mass  distribution:  CH-34,  SA.349/2  and  BO- 105  rotors. 
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1 5  Example  of  CAMRAD  llapwise  moment  error,  CH-34  rotor:  fi  =  0.39,  shaft  angle  =  - 
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f'ig,  16  Components  of  flapwisc  bending  moment  based  on  BRAL  solution,  CH-34  rotor:  /i  =  0.39,  shaft  angle  =  -3“. 


- BRAL  MOMENT 

- El  •  CURVATURE 

- FORCE  INTEGRATION,  100  INTERVALS 

FORCE  INTEGRATION,  300  INTERVALS 


TIME  HISTORY  r/R  =  0.50  STEADY 


I  ig.  17  Comparison  of  flapwisc  bending  moments,  CH-34  rotor:  /i  =  0,39,  shaft  angle  = 
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Fig.  18  Components  of  fiapwisc  shear  force  based  on  BRAL  solution,  CH-34  rotor;  /l  =  0.39,  shaft  angle  =  -5°. 
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Fig.  19  Comparison  of  fiapwisc  shear  forces,  CH-34  rotor;  fi  =  0.39,  shalt  angle  =  -.‘i  '. 
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F'ig.  20  Comparison  of  llapwisc  shear  forces,  CH-34  rotor:  ^  =  0.39 


